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Global  Geoid  and  Gravity  Anomaly  Predictions  Using 
the  Collocation  and  Point  Mass  Techniques 

1.  INTRODUCTION 

The  second-phase  techniques,  point  mass  and  collocation  with  noise,  allow  fine 
structure  of  t he  earth's  gravity  field  to  be  added  to  its  long -wavelength  features 
computed  previously  from  satellite  altimetry.  The  initial,  or  first-phase  treatment 
of  such  data  has  been  carried  out  at  AFGL  in  a  global  short-arc  adjustment  of 
spherical-harmonic  potential  coefficients  (typically  through  the  degree  and  order 
14,  14),  state  vector  parameters  and,  optionally,  certain  tidal  parameters.  This 
adjustment  has  been  the  subject  of  several  AFGL  reports  and  papers,  and  its 
description  will  not  be  repeated  here.  The  most  recent  reference  giving  sufficient 
information  about  its  aspects  is  Reference  1. 

Initially,  a  second-phase  solution  was  conceived  in  terms  of  point-mass 
(P.  M.  )  adjustment  as  the  basic  method  for  a  detailed  resolution  of  the  earth's 
gravity  field  and  its  fundamental  surface,  the  geoid.  Due  to  the  banded  structure 
of  the  matrix  of  normal  equations,  the  P.  M.  parameters  can  be  resolved  in  over¬ 
lapping  strips  eventually  covering  the  whole  oceanic  geoid.  The  P.  M.  approach 
has  been  presented  in  various  AFGL,  reports  and  papers,  most  recently  in 

(Received  for  publication  3  1  August  1984) 

1.  Hlaha,  G.  ,  Hadgigeorge,  G.  ,  and  Rooney,  T.  ( 1984)  Inclusion  of  tidal 
parameters  in  satellite  altimetry  adjustment  model.  Paper  accepted 
for  publication  in  Marine  Geodesy,  Crane  Russak,  New  York. 
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Reference  2  featuring  also  the  results  of  a  global  adjustment.  Its  basic  charac¬ 
teristics  will  be  recapitulated  here  in  a  brief  fashion. 

The  point -mass  approach  has  recently  been  complemented  by  the  method  of 
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least-squares  collocation  with  noise  described  by  Blaha,  which  provides  us  with 
another  technique  to  recover  the  earth's  gravity  field  on  a  local,  regional  or  global 
scale.  In  order  to  achieve  a  detailed  resolution  of  the  gravity  field,  one  should 
proceed  to  an  evaluation  of  a  large-scale  second-phase  solution  in  terms  of  both 
the  point-mass  and  the  collocation  approach.  Accordingly,  this  report  features 
comparisons  between  these  two  techniques. 

As  has  been  indicated  above,  both  second-phase  approaches  are  based  on  the 
first-phase  global  adjustment  of  satellite  altimetry.  In  particular,  the  altimeter 
residuals  from  the  first  phase  are  stored  on  a  magnetic  tape  and  are  utilized  in 
the  second  phase  in  the  role  of  (new)  observations.  Bither  of  the  second-phase 
approaches  accommodate  the  previously  unmodeled  effects  represented  by  these 
residuals.  Thus,  the  high  -  resolution  altimeter  information  enters  the  second- 
phase  adjustment. 

2.  POINT  MASS  APPROACH 

The  P.  M.  adjustment  is  built  on  the  adjusted  geoid  from  the  first  phase  which 

can  represent  a  "normal  field".  The  parameters  in  this  method  are  the  P.  M. 

magnitudes  associated  with  point  masses  distributed  in  an  equilateral  grid.  The 

point  masses  can  form  a  single  laver,  in  which  case  they  are  all  at  the  same  depth 

below  the  surface  of  the  reference  ellipsoid.  Or  the  point  masses  can  form  a 

double  layer  with  twin  point  masses  at  each  location  of  the  grid  separated  by  a 

nredeterm  ined  vertical  distance.  The  number  of  parameters  is  the  same  in  both 

these  modes  since  a  twin  point  mass  represents  only  one  parameter;  the  magnitudes 

of  the  shallower  and  deeper  point  masses  are  equal,  only  their  signs  differ.  The 

resolution  power  depends  primarily  on  the  grid  interval,  and  to  a  much  lesser 

degree  on  their  depth  and/or  vertical  separation  of  the  twin  point  masses.  This 
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is  illustrated  in  1  able  1  of  Blaha  et  al“  presenting  results  from  a  large  test  area 
in  the  southern-most  portion  of  the  Pacific  Ocean.  Thus,  for  example,  if  the  point 
masses  (single  or  twin)  are  distributed  in  a  23  <  2°  equilateral  grid,  the  resolution 
.  orresponds  verv  approximate! v  to  a  <r>0,  00)  spherical-harmonic  expansion. 


2.  Blaha,  G.  ,  Bessette,  R.P.  ,  and  lladgigeorge,  G.  ( 1  !'B4  )  Global  point-mass 
adjustment  of  the  oceanic  geoid  based  on  satellite  altimetrv.  Paper 
submitted  for  publication  in  Marine  Geodesv. 

2.  Blaha,  G.  (ld}34)  First  -  and  Second  -  Phase  G  ravit  v  Field  Solutions  Based  on 
Satellite  Altimetrv.  A  FG  I . -T  R  -  84  -0083  ,  AI)  A14225G. 
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The  second-phase  P.  M.  adjustment  treats  the  new  observations,  the  P.  M. 
parameters  and,  optionally,  certain  tidal  parameters  in  a  simultaneous  least- 
squares  process  leading  to  a  more  detailed  resolution  of  the  earth's  gravity  field 
on  the  global  oceanic  scale. 

The  final  outcome  of  the  two  adjustment  phases  with  respect  to  the  geoid  and 
the  earth's  gravity  field  consists  of  a  set  of  S.  H.  potential  coefficients  and  the 
magnitudes  of  the  point  masses  introduced  at  the  predetermined  locations,  which 
enable  one  to  determine  geoid  undulations,  gravity  anomalies,  and  so  on,  at  any 
points  of  interest.  Rach  such  quantity  is  an  algebraic  sum  of  two  parts,  the  first 
due  to  the  S.  H.  coefficients  and  the  second  due  to  the  P.  M.  parameters.  The 
numerical  values  of  such  quantities  are  considered  to  describe  the  earth's  gravity 
field  to  within  the  desired  resolution. 

3.  COLLOCATION  APPROACH 

The  basic  difference  distinguishing  this  approach  from  the  P.  M.  adjustment, 
from  the  computational  point  of  view,  is  that  a  simultaneous  least -squares  adjust¬ 
ment  of  several  variables  does  not  take  place.  In  principle,  only  one  prediction 
point  is  solved  for  at  a  time.  Therefore,  instead  of  an  inversion  of  a  large, 
strongly  patterned,  system  of  normal  equations,  inversions  of  small  matrices  take 
place.  The  prediction  points  themselves  can  be  distributed  in  an  equilateral  grid 
similar  (or  idential)  to  the  P.  M.  grid.  In  this  case  the  resolution  power,  corres¬ 
ponding  to  an  (n'.n'l  S.  H.  expansion,  is  similar  in  both  approaches. 

The  (new)  observations  are  minus  the  geoidal  residuals  from  the  first  phase, 
and  thev  can  be  limited  to  those  located  within  a  given  spherical  cap  from  the 
desired  prediction  point.  The  predictions  of  greatest  interest  are  those  of  geoid 
undulations,  obtained  through  the  use  of  the  geoidal  covariance  function;  predictions 
>f  other  quantities  such  as  gravity  anomalies  can  be  obtained  for  the  same  location 
upon  using  the  appropriate  cross -covariance  function.  All  the  predicted  quantities 
refer  to  the  "normal  field"  represented  by  an  (n,  n)  S.  H.  expansion. 

Sin  'e  the  mini  tier  of  observations  involved  with  one  prediction  point  can  be 

made  r  easonably  small  (see,  for  example,  weighted  averaging  of  observations 
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locate  |  close  together  as  explained  in  Section  6.  4  of  Rlaha,  )  one  is  faced  with 
inverting  \  small  matrices  as  opposed  to  inverting  an  N'XN  system  of  normal 
equation-;,  where  \  is  the  number  of  prediction  points.  This  allows  for  modest 
.outer  storage  requirements  which  are  set  to  accommodate,  in  one  "loop",  a 
h'iseii  number  of  nre.liet  ion  points  together  with  the  corresponding  vectors  and 
•'  at  rir  es  >t  small  dimensions.  Since  the  prediction  procedure  applies  to  one  point 
at  a  time,  predictions  in  subsequent  loops"  can  he  made  (and  stored  on  a 
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magnetic  tape)  for  the  whole  oceanic  geoid  essentially  in  a  grid  of  an  arbitrarily 
small  interval. 

The  predictions  of  the  (n',n')  geoid  based  on  the  (n,n)  "normal  field"  can  be 
made  in  the  original  equilateral  grid  of  prediction  points,  rather  than  in  a  finer 
grid.  The  original  grid  can  later  be  densified  at  will  upon  using  the  economical 
"errorless  collocation".  Such  a  procedure  utilizes  the  original  predictions  in  the 
role  of  observations. 

The  densified  grid  corresponding  to  the  In',  n')  S.  H.  expansion  can  be  used  for 
a  varietv  of  purposes.  For  example,  it  can  be  used  for  the  construction  of  a 
gooidal  contour  map  (in  this  case  it  may  be  useful  to  make  it  a  fairly  dense  geo¬ 
graphical  grid).  As  another  example,  it  can  serve  in  the  computation  of  the  S.  H. 
potential  coefficients,  via  the  integral  formulas,  through  the  degree  and  order 
(n',n').  A  lower  degree  and  order  expansion,  if  desired,  can  be  obtained  simply 
bv  a  truncation  of  the  expansion  just  formulated.  In  a  very  similar  procedure,  the 
original  as  well  as  the  densified  grid  can  be  made  to  represent  other  geophysical 
quantities  besides  geoid  undulations,  such  as  gravity  anomalies. 


4.  PRACTICAL  APPLICATION  OF  THE  COLLOCATION  ALGORITHM 

We  have  already  pointed  out  the  basic,  or  original,  predictions  are  considered 
to  form  an  equilateral  grid.  Such  a  distribution  yields  a  uniform  resolution  depend- 
ding  on  the  grid  interval  (the  degree  n'  defined  previously  is  assumed  to  be  com¬ 
patible  with  this  interval!.  Accordingly,  a  unit  sphere  representing  the  globe  is 
thought  of  as  subdivided  into  compartments,  each  characterized  by  the  appropriate 
angular  intervals  A<p  andAA,  and  each  associated  with  one  centrally  located  pre¬ 
diction  point.  The  value  of  A<p  is  independent  of  position,  but  the  value  of  AA 
varies  as  A  p/cos  <p,  where  <p  is  the  geocentric  latitude.  The  poles  would  be  an 
exception  to  this  rule  (the  compartments  would  become  spherical  caps),  but  they 
are  omitted  from  this  development.  In  the  case  of  SEASAT,  for  example,  the  ground 
tracks  reach  only  the  latitudes  ±72°  beyond  which  no  predictions  are  made. 

Practical  considerations  coupled  with  numerical  evaluations  have  indicated 
that  including  observation  points  located  within  the  spherical  cap  that  just  covers 
the  corners  of  the  compartment  of  the  prediction  point  called  J  is  reasonable.  This 
provides  for  a  modest  measure  of  continuity  between  neighboring  predictions  because 
of  the*  overlap  in  peripheral  observations.  The  radius  of  such  a  cap,  centered  at  ,1, 
is  approximately  0.  7.7  Ao.  Even  so,  the  number  of  observation  points  involved  in 
one  prediction  may  sometimes  be  prohibitive.  This  occurs  in  spite  of  the  input  data 
being  limited  to  0.  .7  -intervals  along  satellite  arcs,  and  is  due  to  a  substantial 
increase  in  the  number  of  observations  brought  about  by  the  presence  of 
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repeating  tracks.  It  then  becomes  necessary  to  stipulate  a  maximum  allow  ole 
number  of  observation  points  and  discard  the  points  farthest  from  .1  until  this  limit 
is  satisfied,  l  or  this  reason  it  is  useful  to  have  the  observation  points  ordered 
according  to  their  angular  distances  ('i'l  from  .1. 

The  global  collocation  results  are  based  on  the  first-phase  adjustment  of 
S  MAS  AT  altimeter  data  carried  out  in  terms  of  a  (14,  141  spherical-harmonic  model. 
The  original  prediction  points  are  distributed  in  a  2°  X  2°  equilateral  grid  extending 
o'er  the  world's  oceans.  This  corresponds  very  approximately  to  an  S.  11.  expansion 
(UO,  110),  hence  n'  -  !)0.  The  "one  sigma"  of  the  noise  in  the  original  predictions 
has  been  stipulated  to  be  1  m.  This  noise  level  encompasses  the  unaccounted  for 
sea  surface  effects  as  well.  The  densifieation  in  view  of  contour  maps  (for  both 
geoid  undulations  and  gravity  anomalies)  have  been  made  in  a  1°  X  1°  geographical 
grid  using  the  errorless  collocation  technique.  Residuals  have  been  computed  for 
20,  11*2  observation  points  selected  randomly  over  the  global  ocean.  The  root 
mean  square  (rms)  of  the  20,  192  residuals  has  been  computed  as  0.  91  m.  By 
comparison,  the  rms  residual  obtained  through  the  point-mass  technique  has  been 
reported  in  Blah.i  et  al^  as  1.  0  m. 

The  contour  maps  for  N  and  A  g  are  presented  in  Figures  1  and  2,  respectively. 
1'heir  resemblence  to  the  corresponding  contour  maps  in  the  point-mass  approach, 

1  igures  <  and  4,  is  very  satisfactory.  The  quantitative  evaluation  of  this  similarity 
will  be  offered  in  the  next  section.  One  notices  somewhat  irregular  behavior  of 
the  contour  lines  not  only  in  trench  areas,  but  especially  in  areas  containing  a  large 
number  of  islands.  The  latter  problem  is  linked  to  altimeter  errors  which  passed 
undetected  through  the  screening  process.  This  has  been  noticed  mainly  in  the 
West  Pacific,  due  to  a  multitude  of  islands  some  of  which  acted  as  a  reflector  of  the 
radar  pulse  during  data  acquisition.  For  example,  errors  of  several  meters  have 
been  traced  down  in  an  area  south  of  American  Samoa.  With  a  more  elaborate 
screening  process  and,  above  all,  with  an  improved  distinction  between  the  reflective 
surfaces  (sandy  beaches,  and  so  on,  vs  water!  expected  in  conjunction  with  future 
satellites,  this  kind  of  difficulty  will  disappear.  But  with  the  exception  of  the  West 
Pacific,  the  SKASAT  altimeter  measurements  appear  to  be  reliable  over  the  global 
oceanic  surface. 

The  above  discussion  is  illustrated  and  complemented  with  the  aid  of  Table  1 
featuring  the  outcome  of  a  simple  analysis  of  residuals.  These  residuals  account 
for  only  a  small  percentage  of  all  the  available  data,  but  they  have  been  selected 
in  a  random  fashion  to  adequately  represent  the  oceanic  basins.  The  table  is  self- 
explanatory.  fine  notices  a  higher  rms  value  (also  the  average  magnitude)  in  the 
West  Pacific,  the  cause  of  which  has  been  mentioned  in  the  preceding  paragraph. 

By  contrast,  the  Fast  Pacific  exhibits  an  unusually  small  rms  residual.  This  ocean 
b  isin  is  not  only  void  of  major  island  clusters,  but  the  geoidal  relief  is  relatively 
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smooth.  The  W  cst  and  Kast  Pacific  have  been  separated  in  the  computation  by  the 
meridian  of  20a''  longitude. 

If  both  the  Kast  and  West  Pacific.,  representing  the  extreme  values  in  Table  1, 
a  ere  left  out  of  the  computations,  the  results  would  be  (in  the  same  order): 

8,071;  0.05  m  0.03  m,  0.67  m. 

<  >ne  notices  that  the  contribution  due  to  the  Kast  Pacific  lias  had  a  beneficial  effect 
on  all  three  statistics  listed  in  Table  1.  if  no  ocean  basin  should  be  excluded,  but 
tin-  number  of  residuals  in  the  Kast  Pacific  should  be  lowered  to  4,  500  so  as  to 
correspond  more  closely  to  their  number  in  the  West  Pacific,  the  results  would  be 

17,84  1;  0.05  m,  0.01  m,  0.61  m. 

If  the  measurement  errors  in  the  W  est  Pacific  were  eliminated,  the  rms  residual 
\ould  probably  be  brought  below  the  0.9  m  level.  Hut  at  this  point  the  overall 
re.-ult  can  safely  he  presented  as  0.  9  m  rms  residual.  This  is  about  10  percent 
lover  than  the  rms  residual  in  the  point-mass  approach  where,  however,  different 
residua  I  sets  were  used  in  the  computation. 


Table  1.  Results  Pertaining  to  Randomly  Selected  Residuals  of  the 
(  ol location  App roach 


<  Venn 

Niimbi’  r  of 
Residuals 

rms 

Kesidual 

Average 

Residual 

Average 

Magnitude 

North  Atlantic 

1, 956 

0.  89  ill 

o 

o 

3 

0.66  m 

south  A tlantic 

3,  987 

1.01  ni 

0.  04  m 

0.69  m 

1  ndian 

3,  028 

0.  91  m 

0.  00  m 

0.  65  m 

W  .  st  Pacific 

4,  370 

1,  22  m 

-0.  01  m 

0.  75  m 

Kast  Pacific 

6,  85  1 

0.  56  m 

-0.  01  m 

0.  37  m 

Total 

20, 192 

0.  9  1  m 

0.01  m 

0.  59  m 
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5.  COMPARISON  WITH  THE  POINT-MASS  APPROACH 


2 

The  point-mass  approach  and  results  have  been  presented  in  Blaha  et  al  and 
need  not  be  repeated  here.  A  casual  look  at  the  geoidal  and  gravity  anomaly  contour 
maps  may  not  reveal  any  difference  between  the  point-mass  and  the  collocation 
technique,  both  applied  to  the  same  first-phase  adjustment  of  SEASAT  altimetry  in 
terms  of  a  (14,  14)  spherical-harmonic  model.  A  quantitative  evaluation  of  the 
differences  between  these  two  techniques  will  follow.  However,  one  should  first 
realize  a  few  basic  facts. 

Most  importantly,  a  comparison  is  being  made  between  two  completely  differ¬ 
ent  and  independent  approaches,  one  based  on  a  parameter  adjustment  (point-mass 
magnitudes)  and  the  other  based  on  prediction  formulas  using  the  notion  of  co- 
variance  (and  cross-covariance)  function.  The  locations  of  the  original  prediction 
points  in  the  collocation  approach  are  in  general  different  from  the  locations  of 
point  masses  in  the  P.  M.  approach,  although  in  either  case  such  locations  form  an 
approximate  2°  X  2°  equilateral  grid.  However,  the  pertinent  geophysical  values 
are  compared  at  the  level  of  the  densified  grid,  that  is,  the  1°  X  1°  geographical 
grid,  not  at  the  level  of  the  original  prediction  points  or  observation  points.  Thus, 
prediction  errors  due  to  the  densification  increase  the  uncertainty  of  the  values 
actually  compared;  although  this  statement  refers  to  the  collocation  approach,  a 
similar  argument  can  be  made  for  the  point-mass  approach.  The  net  outcome  of 
these  considerations  points  to  a  conservative  nature  of  the  comparisons.  This  is 
even  more  evident  in  the  case  of  gravity  anomalies  (as  opposed  to  geoid  undulations), 
which  can  be  thought  of  as  "doubly  predicted"  quantities,  both  with  regard  to  the 
function  of  the  disturbing  potential  and  to  the  location  of  the  predictions. 

The  blocks  of  comparison  between  the  two  approaches  are  depicted  in  Figure  5. 
These  blocks,  numbered  1  through  13,  cover  all  of  the  world's  oceans,  usually 
more  than  one  block  per  ocean  basin.  The  blocks  have  been  chosen  in  such  a  way 
as  to  avoid  areas  containing  erroneous  data,  especially  parts  of  the  West  Pacific 
as  discussed  earlier,  as  well  as  areas  of  marked  short-wavelength  features.  The 
latter  often  accompany  major  ocean  trenches.  Since  the  comparison  concerns  the 
2° -resolution,  features  corresponding  to  a  1  °  -  resolution  cannot  be  properly  repre¬ 
sented  by  either  method  and  should  indeed  be  left  out  of  consideration.  The  differ¬ 
ences  computed  in  the  1°  x  1°  geographical  grid  are  formed  in  the  sense  the  colloca¬ 
tion  value  minus  the  P.  M.  value.  Figure  5  lists  the  rms  differences  for  geoid 
undulations  (in  meters)  and  gravity  anomalies  (in  milligals)  in  selected  areas. 

Figures  0  through  3  show  the  contours  of  the  differences  in  the  above  functions 
for-  the  blocks  numbered  4  and  8,  which  are  judged  quite  representative  of  the 
overall  comparison. 
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Figure  6.  Contours  of  Differences  of  Geoid  Undulations  (m)  Between  Collocation 
and  Point  Mass  Techniques,  Block  4,  Figure  5 

Finally,  Table  2  summarizes  the  numerical  comparisons  for  all  13  blocks. 

Its  column  structure  resembles  Table  1,  except  that  it  contains  one  additional 
column  featuring  the  range  of  the  pertinent  quantities  (N  andAg)  in  each  block, 
as  estimated  very  approximately  upon  consulting  Figures  1  and  2;  this  range  thus 
contains  the  contribution  of  the  spherical-harmonic  part  as  well.  The  total  rms 
difference  is  0.45  m  for  geoid  undulations  and  2.  8  mgal  for  gravity  anomalies. 
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Figure  7.  Contours  of  Differences  of  Gravity  Anomalies  (mgal)  Between 
Collocation  and  Point  Mass  Techniques,  Block  4,  Figure  5 

The  former  amounts  to  about  50  percent  of  the  one  sigma  in  the  collocation  approach 
alone,  while  the  latter  is  only  about  20  percent  of  the  pertinent  one  sigma. 

If  the  results  of  the  two  approaches  were  random  with  respect  to  egich  other, 
while  each  still  corresponding  to  the  2°  resolution  individually,  such  rms  differ¬ 
ences  would  be  expected  to  be  on  the  order  yj  times  the  appropriate  sigma,  that 
is,  about  1.  2  m  for  N  and  about  18  mgal  for  Ag.  Accordingly,  the  compared  values 
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LONG  I TUDE 

Figure  8.  Contours  of  Differences  of  Geoid  Undulations  (m)  Between  Collocation 
and  Point  Mass  Techniques,  Block  8,  Figure  5 

point  to  an  excellent  agreement  between  the  two  approaches.  In  fact,  the  agree¬ 
ment  with  regard  to  gravity  anomalies  is  better  than  expected,  due  perhaps  to  a 
relatively  narrow  range  of  the  Ag  values  themselves  within  the  blocks  of  comparison. 
One  can  convey  the  idea  of  similarity  between  the  two  approaches  by  stating  that 
the  average  magnitude  of  the  geoidal  difference  is  merely  33  cm  and  the  average 
magnitude  of  the  gravity  anomaly  difference  is  2.  0  mgal,  while  the  average  differ¬ 
ences  are  1  cm  and  0.  08  mgal,  respectively. 
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Figure  9.  Contours  of  Differences  of  Gravity  Anomalies  (mgal)  Between 
Collocation  and  Point  Mass  Techniques,  Block  8,  Figure  5 
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Table  2.  Differences  in  N  and  Ag  Between  the  Two  Approaches 
Compared.  The  upper  values  correspond  to  N  and  are  given  in 
meters;  the  lower  values  (in  parentheses)  correspond  toAg  and 
are  given  in  milligals 


Block 

Number  of 
Points 

rms 

Difference 

Average 

Difference 

Average 

Magnitude 

Range 

1 

455 

0.  56 

-0.  01 

0.  45 

+  20, 

+  65 

(3.  68) 

(-0.  58) 

(2.  86) 

(+  5, 

+  40' 

2 

675 

0.  68 

-0.  02 

0.  52 

-35, 

+  55 

(3.  85' 

(-0.  34) 

(3.  03) 

(-20, 

+  30' 

3 

345 

0.  72 

0.  01 

0.  58 

-35. 

+  25 

(5.  05' 

(-0.  01) 

(3.  57' 

(-20, 

+  20) 

4 

1,410 

0.  43 

0.  01 

0.  33 

-  5. 

+  25 

(2.  52' 

(  0.  10) 

(  1.  97) 

(-10, 

+  10) 

5 

1.131 

0.  56 

0.  03 

0.  43 

-95, 

+  5 

(3.  51) 

(  0.  26' 

(2.  62) 

(-45, 

+  10) 

6 

1,  729 

0.  46 

0.  04 

0.  33 

-35. 

+  45 

.2.  68> 

(  0.  22) 

(  1.  96' 

(-30, 

+  25) 

7 

697 

0.  52 

0.  02 

0.  40 

-15. 

+  5 

13.  61) 

(-0.  26' 

(2.66) 

(-15, 

+  5) 

8 

903 

0.  48 

0.  02 

0.  38 

+  5, 

+  55 

(3.  24) 

(  0.  18) 

(2.  50' 

(-  5, 

+  10) 

9 

437 

0.  33 

-0.  02 

0.  24 

-40. 

-10 

(  1.  80' 

(-0.  17) 

(1.  37) 

(-20, 

+  5) 

10 

615 

0.  32 

-0.  03 

0.  23 

-45, 

+  10 

(  1.  36) 

(-0.  52' 

( 1.  04) 

(-20, 

+  10) 

1 1 

1,  235 

0.  25 

0.  00 

0.  19 

-20, 

+  15 

(  1.  94) 

(  0.  14) 

(1.  39) 

(-10, 

+  10) 

12 

1, 349 

0.  26 

0.  01 

0.  20 

-10, 

+  10 

(1.73) 

1  0.  18) 

(1.  34) 

(-10, 

+  5) 

13 

1,953 

0.  40 

0.  01 

0.  30 

-55. 

0 

(2.  68) 

(  0.  21' 

(2.  03) 

(-20, 

+  10' 

Total 

12.  934 

0.  45 
(2.  84) 

0.  01 
(  0.  08' 

m 

-95, 

(-45. 

+  65 
+  40' 

6.  SUMMARY  AND  CONCLUSIONS 


The  gl  >bal  fi  i  st  -  phase  adjustment  of  satellite  altimetry  described  in  Blaha 
ct  al*  is  a  basis  for  a  more  detailed  second-phase  treatment  of  altimeter  informa¬ 
tion  via  the  point-mass  and  collocation  techniques.  Both  of  the  second-phase  tech¬ 
niques  allow  the  short -wavelength  detail  to  be  added  to  the  long -wavelength  features 
of  tiie  earth's  gravity  field  as  expressed  via  spherical -harmonic  potential  coefficients 
from  the  first-phase  adjustment.  Also  both  approaches  emphasize  the  determina- 
tion  of  the  oceanic  geoid  which  is  directly  related  (except  for  sea  surface  effects) 
to  the  altimeter  measurements.  The  desired  gravity  field  detail  described  in 
this  report  corresponds  to  a  2  °  -  resolution,  approximately  equivalent  to  a  (90,90) 
spherical -harmonic  expansion.  A  brief  summary  of  the  treatment  of  altimeter  data 
including  the  comparison  between  the  two  approaches  follows. 

•  The  point -mass  adjustment  has  proceeded  in  a  single -layer  mode; 
its  main  outcome  (geoid  undulation  and  gravity  anomaly  contour 
maps)  is  presented  in  Figures  3  and  4. 

•  This  mode  is  more  economical  than  the  double -layer  mode,  yet 
the  resolution  power  is  comparable.  ^ 

•The  rms  residual  from  the  point-mass  adjustment  has  been 
reported  as  1.0  m,  which  corresponds  very  well  to  the 
theoretical  sigma  characterizing  the  (90.  90)  truncation 
of  the  spherical -harmonic  mode. 

‘The  smoothed  out  predictions  of  geoid  undulations  and 
gravity  anomalies  in  the  collocation  approach  are  obtained 
in  a  2°  X  2°  equilateral  grid  with  the  aid  of  the  modified 
collocation  with  noise  technique. 

•The  densification  of  this  grid  is  carried  out  by  means  of 
the  errorless  collocation,  resulting  in  a  1°  X  1°  geographical 
grid  suitable  for  the  construction  of  geoid  undulation  and 
gravity  anomaly  contour  maps. 

•The  residuals  in  the  collocation  approach  are  computed 
via  the  errorless  collocation;  this  process  is  compatible 
with  the  construction  of  the  geoidal  contour  map. 

•The  rms  residual  from  a  random  sample  of  20,  192  residuals 
from  collocation  approach  is  computed  as  0.91  m,  which 
accords  very  well  with  the  truncation  of  an  equivalent 
spherical-harmonic  model,  with  the  rms  residual  of  the 
point-mass  approach,  and  with  the  one  sigma  associated 
w  ith  the  geoidal  contour  map. 

‘The  average  residual  is  similarly  computed  as  1  cm  and 
the  average  magnitude  of  the  residuals  is  computed  as 
59  cm. 

•  The  comparison  of  results  for  both  and  Ag  values  proceeds 
in  13  blocks  depicted  in  Figure  5,  encompassing  12,934 

grid  points. 
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•The  mis  .lilfereme  in  N  between  the  two  techniques  is 
computed  as  0.45  m;  tlie  rms  difference  in  Ag  is  2.8  mgal. 

•The  average  difference  in  \  is  !  cm.  while  it  is 
0.08  mgal  for  A  g:  the  average  magnitudes  at' the  differences 
are  22  cm  and  2.0  mgal.  respectively. 

•  The  geoidal  and  gravity  anomaly  results  form  the  two 
approaches  are  thus  seen  to  correspond  very  closely  to 
each  other.  It  can  be  concluded  that  these  techniques 
nnl'irm  and  verify  each  other. 


END 

FILMED 

7-85 

DTIC 


